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Abstract: We compared changes in the spatial clustering of schistosomiasis in Southwest 
China at the conclusion of and six years following the end of the World Bank Loan Project 
(WBLP), the control strategy of which was focused on the large-scale use of chemotherapy. 
Parasitological data were obtained through standardized surveys conducted in 1999-2001 
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and again in 2007-2008. Two alternate spatial cluster methods were used to identify spatial 
clusters of cases: Anselin's Local Moran's / test and Kulldorff s spatial scan statistic. 
Substantial reductions in the burden of schistosomiasis were found after the end of the 
WBLP, but the spatial extent of schistosomiasis was not reduced across the study area. 
Spatial clusters continued to occur in three regions: Chengdu Plain, Yangtze River Valley, 
and Lancang River Valley during the two periods, and regularly involved five counties. 
These findings suggest that despite impressive reductions in burden, the hilly and 
mountainous regions of Southwest China remain at risk of schistosome re-emergence. Our 
results help to highlight specific locations where integrated control programs can focus to 
speed the elimination of schistosomiasis in China. 

Keywords: schistosomiasis; spatial pattern; clustering; the World Bank Loan Project 
(WBLP); hilly and mountainous regions; Southwest China 



1. Introduction 

Schistosomiasis remains a major public health problem in many developing countries in tropical and 
subtropical regions [1], with more than 200 million people infected and 779 million people at risk [2]. 
The global burden of schistosomiasis has been estimated at 1.7 to 4.5 million disability-adjusted 
life years (DALYs) [2], but the true burden is proved to be between four to 30 times greater than 
previously expected [3]. Of the three main schistosome species, Schistosoma japonicum is responsible 
for human and animal infections in southern China, the Philippines, and parts of Indonesia [4]. 

China was one of first countries in East Asia to initiate a national schistosomiasis control program. 
Substantial achievements in fighting against this disease have been made since the mid-1950s [5]. 
By the end of the 1980s, schistosomiasis was eliminated in four provinces, but the disease remained 
endemic in eight others [6,7]. In a renewed attempt to eradicate schistosomiasis in China, a 10-year 
control project was launched in 1992 with a World Bank loan of US$ 71 million and a supplementary 
fund of US$ 82 million from the Chinese government [8]. This aggressive project employed large-scale 
chemotherapy, and also included additional interventions such as health education, chemical control of 
snails and other environmental exposure modifications to reduce exposure. An evaluation of this World 
Bank Loan Project (WBLP) in 2002 [9] found that infection rates in humans and livestock had 
decreased by 55% and 50% respectively and that the number of acute infections and of individuals 
with advanced disease had also significantly decreased. One province that took part in the project, 
Zhejiang Province (one of the remaining eight endemic provinces), was able to fulfill the national 
criteria for elimination. 

Although clear progress in schistosomiasis control was made during the period of WBLP period, the 
long-term progress toward elimination has been less clear. For example, in recent years there have 
been reports of the reemergence of schistosomiasis in some intervention areas [10]. Potential 
contributors to schistosomiasis reemergence include environmental and socio-political changes, allowing 
for the establishment of new snail habitats and subsequent rebounds in disease burden [11,12]. 
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Given concern about the durability of the progress against schistosomiasis achieved during the WBLP 
period, we sought to assess the longer-term progress after the end of this project. 

In this study, we investigate temporal changes in the spatial clustering patterns of schistosomiasis in 
Southwest China to better understand the long-term progress toward eradication after the end of the 
WBLP. Data from the national schistosomiasis disease registry are used to map the prevalence 
of schistosomiasis during two periods of time for which these data were available 1999-2001 and 
2007-2008. We focus our analyses within the two hilly and mountainous provinces in China in which 
5*. japonicum was endemic and where the WBLP interventions were implemented. The first period 
(1999-2001) allows us to identify spatial disease patterns near the conclusion of the WBLP, while the 
second period (2007-2008) permits investigation of longer-term changes in disease clustering after the 
end of the project. 

2. Materials and Methods 

2.1. Approach and Study Area 

Techniques of geographic information system (GIS) and methods of cluster detection were 
combined to evaluate temporal changes in the spatial distribution of S. japonicum. The analysis was 
conducted at the county level within the two mountainous provinces of Southwest China, Sichuan and 
Yunnan where the WBLP was focused (Figure 1). Sichuan ranks as both the 3rd largest (485,000 
square kilometers) and 4th most populous province (80.4 million in 2010) in China. Adjacent to 
Sichuan, Yunnan is a province spanning approximately 394,000 square kilometers and with a 
population of 45.9 million (2010). Most of the both provinces are mountainous or hilly, with plains and 
plateaus covering the rest of the land area. 

Figure 1. Location of Sichuan and Yunnan Province. Geographical rivers layer was 
overlaid. The map was created using the ArcGIS software (version 10.0, ESPJ Inc., 
Redlands, CA, USA). 
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2.2. Parasitological Data 

National schistosomiasis datasets for the two analysis periods (1999-2001 and 2007-2008) were 
provided by Fudan University (formerly Shanghai Medical University) and the National Institute of 
Parasitic Diseases in Shanghai, Chinese Center for Disease Control and Prevention. The databases are 
county-based, with all reported schistosomiasis cases and the population at risk in each county 
included. Additional details describing the process of data collection can be found elsewhere [13]. 

2.3. Statistical Analysis 

Preliminary description of prevalence of schistosomiasis. Mapping of prevalence of schistosomiasis 
for each year was performed to assess the spatial distribution of the epidemic over time. Spatial 
autocorrelation was investigated using the Global Moran's / statistic [14]. 

Detection of spatial clusters. For each year, spatial clustering was assessed using Local Moran's / [15] 
and Kulldorff statistic [16] to identify high-risk areas. These two methods provide complementary 
information on clustering patterns to identify areas of relatively high prevalence. The Local Moran's / 
statistic provided a measure of the spatial autocorrelation for a given county with its neighbors. For 
both Global and Local Moran's /, a pair of neighbor counties is defined with those sharing a common 
border. We assumed the weight element as 1 if two counties are neighbor and 0 otherwise in 
the weight matrix which defines the spatial relationships among counties. Spatial clusters of 
schistosomiasis cases are identified by finding local areas where high prevalence counties border other 
high prevalence counties (defined as high-high county) and where high counties border low prevalence 
counties (defined as high-low county). The Kulldorff s statistic uses a moving window of pre-specified 
shape to identify single counties or groups of counties with statistically significant relatively high risk. 
We used an eclipse window in this study which moved throughout the study area. 

Frequency of occurrence of spatial clusters. The frequency of spatial cluster occurrence was 
defined as the number of years during which the county contributed to a cluster and this was calculated 
for each county during the two study periods. Two types of frequencies were considered: (i) clusters 
detected by one of the two statistics above (weak evidence of clustering) and (ii) clusters detected by 
both statistics above (strong evidence of clustering). Mapping of prevalence of schistosomiasis was 
performed using ArcGIS software (version 10.0, ESPJ Inc. Redlands, CA, USA) and the Global and 
Local Moran's / statistic were calculated using the ArcGIS software as well. Spatial cluster analyses 
with Kulldorff s statistic were implemented using the SaTScan software (version 8.0, Kulldorff and 
Information Management Services, Inc., Boston, MA, USA). 

3. Results 

Figure 2 shows the annual prevalence of schistosomiasis during the two study periods, from which 
we can see that the burden of schistosomiasis decreased dramatically from the first (1999-2001) to the 
second (2007-2008) period. Actually, the mean prevalence of schistosomiasis decreased dramatically 
from 100.5 cases per 10,000 people in 1999 to 6.02/10,000 in 2008 during the two study periods 
(Table 1). This downward trend was accompanied by decreasing variation in prevalence across counties 
with the interquartile range (IQR) contracting from 0-162.2/10,000 in 1999 to 0-7.63/10,000 in 2008. 
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The Global Moran's / coefficient was particularly low with p > 0.05 for each year except 1999, 
indicating no significant overall spatial autocorrelation. 

Figure 2. Annual prevalence of schistosomiasis at county level during the two study periods. 
1999 2007 




The maps were created using the ArcGIS software (version 10.0, ESRI Inc., Redlands, CA, USA). 
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Table 1. Prevalence (cases/10,000) of schistosomiasis and its global spatial autocorrelation 
in the mountainous area, China. 



Year 


Min 


Qi 


Mean 


Q3 


Max 


J 


P 


1999 


0 


0 


100.5 


162.2 


1172 


0.13 


0.04 


2000 


0 


0 


46.98 


10.78 


672 


-0.03 


0.55 


2001 


0 


0 


51.66 


36.22 


595.5 


0.07 


0.16 


2007 


0 


0 


6.85 


5.69 


70.73 


0.11 


0.07 


2008 


0 


0 


6.02 


7.63 


29.84 


0.09 


0.13 



/: the Global Moran's / coefficient; p: /j-value for the Global Moran's / statistic. 



Figure 3. Annual spatial clusters of schistosomiasis during the two study periods. 




Each panel shows for each year the results of both methods, the Anselin's Local 
Moran's / test and the Kulldorff 's spatial scan statistic. The maps were created 
using the ArcGIS software (version 10.0, ESRI Inc., Redlands, CA, USA). 
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Spatial clusters are shown in Figure 3, which were detected by both the Kulldorff s spatial scan statistic 
and the Local Moran's / results. The Local Moran's / identified 14 high-high counties during 1999-2001, 
with an annual number ranging from 8 (in 2000 and 2001) to 10 (in 1999) and the median annual 
prevalence of counties flagged as being in high prevalence clusters of 216.79 (IQR = 160.03-450.81) cases 
per 10,000. In 2007-2008, the Local Moran's / identified 18 high-high counties, with an annual 
number ranging from 10 (in 2007) to 17 (in 2008) and the median annual prevalence within these 
counties of 19.52 (IQR = 12.96-27.01) cases per 10,000. Note that no high-low counties were found 
during the two periods. The Kuldorff spatial scan method identified 12 different significant spatial 
clusters during 1999-2001 (four in each year); the median number of counties included in each cluster 
in this early period was 3 (IQR = 1-10). In 2007-2008, the Kuldorff spatial scan statistic identified 
13 significant spatial clusters altogether (5 in 2007 and 8 in 2008); in this latter period the median 
number of counties contained within a cluster dropped to 1 (IQR = 1-3). Over the two study periods, 
almost 88% of the high-high counties identified with the Local Moran's / were included in clusters 
detected by SaTScan and 64% of the SaTScan clusters encompassed high-high counties. 

Spatial clusters generally occurred mainly in three regions (Chengdu Plain, Yangtze River Valley, 
and Lancang River Valley), but in different counties from year to year over the two study periods 
(Figure 4). Among the counties contributing to a cluster at least once, the median frequency was 2 
(range = 1-5) for clusters detected by at least one method, and 1 (range = 1^1) for clusters detected by 
both methods. Five counties were detected three or more times by both methods (Renshou, Meishan, 
Dechang, Eryuan, and Dali). 



Figure 4. Frequency of cluster occurrence during the two study periods. 
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(A) Frequencies of occurrence of spatial clusters detected by one of the two methods (Anselin's 
Local Moran's / test or Kulldorff 's spatial scan statistic); (B) Frequencies of occurrence of spatial 
clusters detected by both methods. The maps were created using the ArcGIS software (version 10.0, 
ESRI Inc., Redlands, CA, USA). 
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4. Discussion 

This study presents an application of spatial cluster detection methods to examine changing patterns 
in the burden of schistosomiasis in hilly and mountainous areas in Southwest China during the period 
immediately following and 6-7 years after the end of the WBLP. Previous studies [8,17,18] have 
reported the evaluation of impact of the WBLP control strategies on schistosomiasis control, but most 
of these studies evaluated only the short-term impact of the WBLP, and did so from a non-spatial 
perspective [19]. Since different patterns of spatial distribution would require distinct control 
strategies, our study provides new insights and valuable information for future planning of 
schistosomiasis control efforts. To our knowledge, only Zhang et al. [13] have so far investigated 
spatial variations of schistosomiasis risk in this area using a Bayesian hierarchical model, but this 
analysis did not attempt to identify areas of spatial clustering. Spatial cluster analysis with time 
sequence provides unique insights since it provides valuable clues about potential disease sources 
(and determinants) [20], areas where clusters may recur [21], and promising areas to target for 
increased control efforts [22,23]. 

It was found that the overall prevalence of schistosomiasis had substantially reduced to a low level 
in the period of 2007-2008. While causality cannot be unambiguously attributed to the WBLP 
activities, we believe it is reasonable to conclude that these efforts may have been responsible for such 
impressive reductions in disease burden. Our analysis showed that there was a significant spatial 
autocorrelation of prevalence of schistosomiasis (I = 0.13, p = 0.04) in 1999 but no significant spatial 
autocorrelation existed in the following years, suggesting that the spatial distribution of the disease 
changed from a clustered pattern to a highly fragmented one. 

Our spatial cluster analysis also demonstrated the greatly reduced burden of disease between these 
two periods: the median annual prevalence of counties within detected clusters decreased from 216.79 
to 19.52 cases per 10,000 (Local Moran's / results) and at a smaller scale (the median number of 
counties per cluster reduced from 3 to 1, Kulldorff s spatial scan statistics). This reduction is evident in 
Figure 2. However, we found in the latter evaluation period that more clusters of high-risk counties 
were widely scattered than in 1999-2001 (Figure 3), indicating there were still large population at risk. 
This suggests that the WBLP strategy may have mitigated the burden of schistosomiasis but failed to 
effectively compress the spatial extent of the endemic areas. 

It is well-known that the transmission of S. japonicum is closely related to the distribution of the 
snail intermediate host, which is associated with environmental conditions such as vegetation coverage, 
land-use patterns, quality and humidity of the soil [24]. In mountainous regions of Southwest China, 
snails are distributed in the paddy fields, irrigation ditches, and their adjacent meadows with water, 
where human and livestock (especially cattle) are frequently found [25]. The scattered clusters of 
schistosomiasis reflected the dispersed distribution of snail habitats to some extent. Possible causes of 
this pattern are environmental management [7], which included replacing farmland with forest or grass 
and turning rice paddies into cotton or corn land in recent years, and development of roads or 
highways with the social development in China [26], some of which passed through the snail habitats 
creating more snail habitats at a smaller scale. 

In general, clusters occurred in different counties from year to year, but were more frequently found 
in five counties. Map A (Figure 4) shows that clusters mainly occurred in three regions: Chengdu Plain, 
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Yangtze River Valley, and Lancang River Valley, indicating these regions should remain the focus for 
schistosomiasis control. Furthermore, the strong evidence of clustering, shown in Map B (Figure 4), 
implies that the five counties (four in yellow and one in orange) should be given priority in future 
schistosomiasis control efforts. 

The control of S. japonica requires interruption of the parasite lifecycle by eliminating the 
intermediate host and the parasite from the definitive host and preventing infection of the intermediate 
or definitive host [27,28]. The provision of praziquantel in highly endemic regions is effective at 
reducing the infection rate [29,30], however, drug treatment alone may not be able to mitigate 
schistosome re-emergence or the establishment of new endemic areas [31]. This is exemplified by 
our study where, despite a decrease in prevalence, clusters became more scattered in 2007-2008. 
Therefore, an integrated control program [31] is required to combat schistosomiasis to achieve the aim 
of reducing human infection and prevalence to below 1% by 2015 in all endemic areas [32]. These 
interventions should include health education, access to clean water and adequate sanitation, 
mechanization of agriculture, fencing of water buffaloes, snail control, chemotherapy, and even future 
vaccination of livestock [6,33,34]. 

There are several limitations to the analysis presented here. First, our analysis depends on the 
reliability of reported surveillance data and stability in the efforts to detect cases over time. We do not 
believe that there have been substantial changes in reporting over the period of this analysis; Second, 
we are unable to unambiguously attribute changes in schistosomiasis prevalence to the WBLP since 
this was not a cluster randomized study and we do not have control areas with which to compare 
spatial patterns of disease burden over time; Third, methods of cluster detection we have used have 
limitations. For instance, the Local Moran's / is biased when population at risk in each county varies 
substantially (on which prevalence is calculated) over the study area [35] and the Kulldorff s spatial 
scan statistic requires that we pre-specify the shape of the clusters to be detected. We combined these 
two different approaches for cluster detection methods, as has been recommended by others [36], to 
help address these limitations and to improve our overall understanding of the complex and changing 
spatial distribution of schistosomiasis. Furthermore, the two methods show the trend of schistosomiasis 
from different aspects. The Local Moran's /, based on spatial attribute (e.g., prevalence), reflects 
severity of schistosomiasis; while the kulldorff s statistics, detecting spatial distribution, focuses on 
extent of schistosomiasis. 

5. Conclusions 

In summary, the period following the WBLP intervention was associated with substantial 
reductions in the burden of schistosomiasis, but limited impact on the spatial distribution of clustering 
within the study area. Clusters mainly occurred in three regions: Chengdu Plain, Yangtze River Valley, 
and Lancang River Valley, and that clusters were most frequently found in five counties, which must 
be targeted for increased interventions to speed the elimination of schistosomiasis. 
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